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Abstract 

Approximate solutions to functional evolution equations are constructed through a combination of series and con- 
jugation methods, and relative errors are estimated. The methods are illustrated, both analytically and numerically, 
by construction of approximate continuous functional iterates for sin a;, and Xx{l — x). Simple functional 

conjugation by these functions, and their inverses, substantially improves the numerical accuracy of formal series 
approximations for their continuous iterates. 




Functional iterates of the sine 
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I. OVERVIEW 



There are several circumstances where one encounters the functional evolution equation [l|, [l3| , 

Xs+t ^XsOXt , (1) 

corresponding to an abelian flow in the subscript parameter, with initial condition Xt=o — id, the identity 
map. Here, "o" is functional composition. For example, such an equation governs the changes in scale 
that underlie the renormalization group or the rectification of local trajectories in dynamical systems 

0. While © does give rise to a local differential equation, in principle, the actual problem of interest may 
instead present the result of a fixed "large" step in t, say xi, from which all continuous iterates Xt are to be 
determined. 

In this situation, a useful direct approach is to construct an A^th order formal series approximation for xt 
around a fixed point of xi, say at x = 0, by series analysis of the s = 1 case of ([T]), written as 

xi {xt (x)) = Xt {xi (x)) . (2) 

With initial conditions corresponding to Xt=o — id, the result has the form 



(A^ approx 

\-^J — 

fc=l 



(a:)=^c,(Ox\ (3) 



where ci (0) = 1, Cfc>i (0) — 0, so x[^q^^^'^°^'' (x) = x. Note that it is only necessary to construct accurate 
approximations to Xt for any unit interval in t, for then the composition rules xi^t — xi o xt and = 
X-i o Xt can be used to reach higher or lower values of t. 

Now, for a given N the series ([3|) may not produce accurate results for values of t, or for an interval of 
initial x, of interest. Quite possibly this may be overcome if the series is convergent, by taking larger N ^ 
but if the series is only a formal one, perhaps asymptotic, with zero radius of convergence, larger N may not 
be a viable option. So what can be done then? 

Formally, for any fixed N , with a;„ = xio---oxi,it follows that 

n compositions 

r {I^ approx) / . 

Xt = lim XnO Xt O X-n (4) 

n— f ±00 

is a solution to This is a consequence of 

/ (N approx) \ / [N approx) \ /r-x 

Xio\XnOXl ' o a;_„ = I a;„+i o ' o a;_„_i o xi . (5) 



In the limit n — > ±00 Eqn ^ is obtained, ij either of those limits exists. For a specified class of problems 
it might be possible to give an elegant proof that either limit exists by using various fixed point theorems 
from functional analysis 0, but that is not our objective here. 

The objective here is to estimate the numerical accuracy obtained by conjugating ([3]) a finite number of 
times, n, with the given, exactly known, "finite step" functions x±\. That is, our concern here is the relative 
error obtained as a function of n, prior to taking the limit n — 00. In many cases, n-fold conjugation with 
the given xi, and its inverse X-i, dramatically improves the numerical accuracy of the series approximation, 
with the error vanishing exponentially in either TV Inn or Nn. This behavior came to light in follow-up 
studies of earlier work 



II. A RATIONAL ILLUSTRATION 

As a tractable example d, consider 



2 



In this particular simple case an exact solution to ([T]) is 

= (7) 

But to illustrate the series method, consider ^ given only xi in Recursive analysis in powers of x 

immediately gives Ck {t) = t^^^ , where we have defined the form and scale of the t parameterization by the 
choice C2 {t) = t, to be in agreement with ([7]). For instance, with Xt {x) = x + tx"^ + C3 {t) x^ + O (a;"*) we 
find 

^""^ ( = - C3 it)) x^ + {x") . (8) 



I - Xt [x) 



Thus (21) is satisfied if and only if C3 {t) — . So it goes to higher orders in x, with each Ck+i determined 
by Cj<k to satisfy ^ — not just by direct expansion of ([7]). 
The result for the approximation is therefore 

^{N approx) = ^ ^fc-l ^fe = ^ t ^X ) X 

k=\ 

Then, by n-fold conjugation of this approximation with x\ and its inverse, obtain 

...o4— o. (l-(^)'')/ + ^ (10) 

This indeed converges to the exact result, ([7]), as ri cxd for any fixed N > 1. But what is the relative error 
for finite n? 

Since we know the exact answer in this elementary case, this error is not difficult to compute. For any 
fixed N and finite n the relative error is 

dcfn. a;t (x) - a;„ o ''■'P™''-' o x_„ (x) 1 - tx + nx f tx \^ 

y 1+nx 1 

For > 1 this indeed vanishes as n 00, for any fixed x, so long as tx ^ 1. However, more importantly, 
for n finite but large compared to both t and 1/x, this error goes to zero as the [N — l)st power of n, 

t^x 1 

Rt (x,N,n) « ^ . (12) 

^ ' ' ' «»|tMl/:i;| l-tx ^ ' 

Therefore, in principle, one can obtain numerical results to any desired accuracy by repeated conjugation of 
the approximate series with the given step- function xi. It is remarkable that it is only required to take any 
fixed TV > 2 to produce such results, although in practice, as is manifest in (fT2|) . computational efficiencies 
can be improved by taking larger N since the desired accuracy is then reached for smaller n. 

The result for the relative error can be understood in terms of Schroder theory [l^l as it applies to this 
simple example. The continuous iterate xt given in ([7]) may be constructed from a Schroder function, ^, 
and its inverse. An exact result for this particular example is given by ;17l] 

=exp(-l/a;) , (a;) = -1/ In (x) , (13) 

In Schroder's conjugacy framework, the general iterate is xt (x) = ^'^^ (e*^ i^)): or, 

^'(xt (x)) = e**(x) . (14) 

Thus w = 'if (x) is just the change of variable that reduces the effect of evolution in t to nothing but a 
multiplicative rescaling, wt = e*w. 

But suppose that the exact xt were supplanted by an approximation of the form x[^ approx) _ xt+ax^^^ + 
O (xf^^^) , for some coefficient a. Then 

* (xf "PP''™^) = (1 + axf -1 + O (xf )) exp (-l/x*) . (15) 
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Alternatively, with xt = — l/ln^(a;() = —l/\nwt, 

.(appro.) ^ ^ UN appro.) + Q ( / ] ] . (16) 

^ ^ V (-Inwt) \(-lni«t) // 



Then it follows from the multiplicative rescaling behavior of w that under the variable change Xt — > Wt the 
conjugated approximation x„ o a;^^'''^™ o a;_„ presents itself as 

W„ o Wt^PP'™ o w_„ = I 1 + , " ,jv_i + O [ ^ ) I ^* • (1'^) 

\ (n-lnwt) Un-lnwf) // 



This gives a relative error with the same power law asymptotics [l^ as ([T2|) . namely, i? ^ a/n 



JV-l 

n—foo 



III. THE ROOTS OF SIN 

A more interesting example is provided by the sine function, where our notation for the continuous iterate 
is now siut (x) with sini (x) = sin (x) and sin_i (x) = arcsin(x). In this notation, the abelian functional 
composition equation is 

siug o sint = siug+t . (18) 

Again specializing to s = 1, written as 

sin (siuf (x)) — sint (sin (x)) , (19) 

we find a formal series solution, 

sini—) (X) = x - itx3 + i^x^ - (iI^!!_^^±M!,7 ^ i2^t-ni^-nft^, ^ ^ (^n) . 

3. 5. 3x7. 9. ^^^^ 

This is sufficient for the numerical work to follow (other results for this example are available online: 
ISchroeder.html' & 'TheRootsOfSin.pdf). 

In this case, the series ([20]) is probably not convergent for almost all t [1, [l^ [l^ , although it obviously is 
convergent for < G Z. Rather, for generic t the series appears to be asymptotic. For example, for t = 1/2, 
using Mathematica to extend the series ([20| to O (x^^) or so, the smooth behavior of the series coefficients 

Cfc for fc < 61 suggests a finite radius of convergence « 4/3, as estimated by 1/ Icfcl^^'^. But then less smooth 
behavior sets in for the Cfc>6i and numerical estimates of the radius begin to fall towards zero, as would be 
expected for an asymptotic series. 

Also note in passing that (|20l) implies the corresponding Schroder function has an essential singularity at 
X = 0, as given explicitly by 17] 

*W^5„.«'=exp(^ + ^.' + J^.. + 0M) , ,21) 

to be compared to (|13p . This follows from ^I^ (x) = exp (/ dx/u (x)), using u (x) = dsuYt{x) /di|f^Q ~ 

fapprox), . . 

t=0 

Nonetheless, the conjugation method produces approximations siut « An^t which provide compelling 
numerical evidence for the existence of a limit, hence an exact result for sint, as n — > oo. But note it is 
important in this case to take the conjugations to be of the form 

■ • ■ (approx) . /r,r,\ 

An^t = sm_„ o smt o sm„ , (22) 

the general rule being to act first with functions that are smaller in magnitude than the identity map, 
thereby leading to evaluation of the truncated series at points closer to x = 0. 
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The improvement wrought by conjugation of the truncated ninth-order series in (|20l) is easily seen in 
the following graphs, for the case t = 1/2. The series itself is not credible beyond x = tt/2, but a single 
conjugation forces the expected periodicity in x, and gives the more plausible green curve seen in the first 
Figure. Repeated conjugation does not produce any discernible differences with the green curve to the 
accuracy of that first Figure, but when the graph is magnified, as it is in the second Figure, numerical 
evidence for convergence of the sequence of conjugations is quite compelling. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 



1.580 1.585 1.590 
X 



The ninth order series ^ for t = 1/2, in black, along The i = 1/2, n = 1, 2, 3, 4, & 5-fold conjugations (|22l) . 
with its n = 1 conjugation ((22|) . in green. in green, orange, blue, red, and purple, respectively. 



Proceeding in this way leads to the set of curves for various values of t shown in the following Figure. 
Each of these curves results from the five-fold conjugation A^ t of the truncated ninth-order series (|20p . 
Note once again, as previously remarked in a general context, it is only necessary to construct accurate 
approximations to siut for any unit interval in t, for then the composition rules sini+t = sinosiuf and 
sin_i_|_f = arcsino siuf can be used to construct the curves for higher or lower values of t. Also note from the 
numerics the obvious inference that sint (x) becomes the periodic triangular "sawtooth" function as t — > 0, 
with limt^o sint ((2fc + 1) 7r/2) = (-1)'' tt/2. 




Various sint {x) as given by five-fold conjugations Ag^t (x) of the ninth-order series ([20|) . for t < 1. 
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As constructed, sin^ (a;) for t > is guaranteed to be real for all real values of x, but certainly it is not 
obvious for generic t what numerical values are actually attained at the extrema for x — iTrmodTr. It 
suffices here for us to point out that the maxima are approximated by the simple expression: 



sint (7r/2) 



(23) 



At least, this is true for < t < 1, where the relative error between the exact (numerical) value of sint (7r/2) 
and this approximation is less than about 3 parts per mille. (More accurate numerics are available online ) 
The branch point at ^ = exhibited in this approximate expression is perhaps the most direct numerical 
evidence that the iterates are not analytic a.i t — for all x. 

The graphs in the second Figure above give a sense of the overall relative error, but lacking closed-form 
expressions for either the iterates or the conjugations of the series approximations, closed-form results for 
the error are not available for generic t. For i G Z, however, precise calculation is indeed possible since 
both exact results and convergent series are known. It suffices here to consider just one exact case, t = \. 
Defining the relative error as before. 



i?i (x, n) 



■ (approx) . / N 

sm X — sin_„ o smj^ o sm„ (x) 



smx 



(24) 



we have computed numerically the error involved in various conjugations of the ninth-order series. As 
previously remarked, conjugation by the sine guarantees the approximations are periodic. The results are 
shown here. 



R 



O.Oe+O 



-l.Oe-9 



-2.0e 




Rt=i {x, n) for n = 4, 5, and 6, in green, red, and blue, respectively. 

It would be interesting to compute relative errors for other, generic but at this stage it is only possible 
for us to compute the relative successive differences^ 



St {x,n) 



. (approx) . / N . . (approx) . / 

sm_„ o smj o sm„ [x) — sm_„-|_i o siUj o sin„_i [x) 



sm_„ o sm. 



(approx) 



For example, consider t = 1/2 using the ninth-order series. 



o sin„ (x) 



(25) 



l.Oe-5 



5.0e-6 



O.Oe+0 




St=i/2 (2^7 "•) foi' n = A,b, and 6, in green, red, and blue, respectively. 

This numerical data supports the proposition that these continuous iterates of the sine function are well- 
defined and straightforward to compute using the methods in this paper. 
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IV. THE LOGISTIC MAP 



The sequence of conjugations converges more rapidly in situations where the underlying Schroder function 
is analytic about the fixed point, in contrast to (fT3|) and (f2T|) . Instead of power law behavior, for such 
situations the relative error vanishes exponentially in n, the number of conjugations. The general theory is 
well-illustrated by the logistic map, 

xi{x) ^ Xx{l - x) , for 0<a;<A/4 and < A < 4. (26) 

The result of the theory is the following: 

[Theorem] Relative error after n-fold conjugation of the truncated series ([3]) is given by 

Xt (x) - Xn {x[^ appiox) ^^^^^^j 

Rt{x,X,N,n) = ^ ^ 



A 



Xt (x) 
r {x, t, A, N, n) 



(27) 



[Proof] To sketch a proof, and to see more clearly the assumptions involved, as well as to obtain expressions 
for r (x, t, A, n), write the truncated series as 

^{N apx) ^ _ x^'+^Sm {x, t) , (28) 

where Sn ix,t) represents the exact difference, whose expansion in x begins 0{x°). Thus the conjugation 
gives exactly 

Xn (x[^ '^^'^^ (X_„ (x))^ = Xn (^Xt (x_„ (x)) - (x_„ (x))"^^ Sn {X-n (x) , t)j . (29) 

Now expand the RHS in powers of (a:_„)^^^ Sn, 

Xn (xt {x^n (x)) ~ {x^n (2^))^^^ (a^-n (a;) , t)^ = x„ {xt (x_„ (a;))) (30) 

- (X_„ [x])^^^ Sn {X-n (x) , t) x'„ {xt {x^n {x))) + O {[x^n)'^^^ S%^ 

Since it consists of exact trajectories that obey ([T]), the first term gives x„ {xt {x-n {x))) = Xt (x), while the 
second term involves 

x'nixtix^nix))) = —1 -^Xn {Xf (X-n (x))) = -^Xt (x) , (31) 

— Xt\X^n{x)) —Xt-n\X) 

ax ax 

again using ([T]). Thus 

Xt{x)-x[''^^^^-^Hx) = {X-n{x)f+' SN{x,n{x),t) ^1^1^ + O {{xlnf S%) . (32) 

To proceed, we require that the unit step function is such that x-n (x) flows toward the fixed point at the 
origin for the point x under consideration (so we suppose |A| > 1, but if not, just interchange xi and X-i). 
We also suppose that n has been chosen large compared to t so that 

xt_„(x) = A*-"£„(x,t) (33) 

is small, where e„ {x,t) = x + O (x'^) . If these requirements are met, then 



Xt {X) - -'^ N) ^ ,-.(.+1)^^+1 0) Sn {X-Sn {X, 0) , t) ^Jn^^^l^'y,, 



O 



A £„ On 



(34) 
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where dSn {x,t) /dx = 1 + O (x). The result for the relative error is therefore of the form in the statement 
of the Theorem, with 

r [x, t, A, N, n) = g» O^a^ (A-"£n (x, 0) , dxt {x) /dx ^ ^ _„(^+2)^2^+2^2^1 _ (35) 

A dSji yx^ t) / dx Xt \Xj L J 

For small x we have 

e^+^ {x, 0) = x^+i (1 + (x)) , den (x, t)ldx^l + {x) , 
Sn (A-"e„ (x, 0) , t) = CN+i {t, X) + (x) , {'^'' ^ - [1 + O (x)) , (36) 

Xt (X) X 

and therefore 

r {x, t, A, N, n) = X'^cn+i {t, A) x^ + O (x^^^) , (37) 

again as previously stated. ■ 

The form given in psp enables analysis of the size of the error as a function of x. Moreover, the form in 
P7p suggests an approximate scaling law for the error 

Rt (x, A, N, n) « A^i?t [x, X,N,n + l) , (38) 

at least for x near the origin. It is interesting to check whether this is true for larger x. In fact, it is. 
To be specific, consider the case A = 2 which can be solved in closed-form to obtain: 

* (x) = -i In (1 - 2x) , (x) = i (1 - , (39) 

(x) = ^ (1 - (1 - 2xf) , ^ = 2* (1 - 2x)-'^'' , (40) 
e„ {X, t) = 2—1 (1 - (1 - 2xf-^) , = (1 - 2xr'^''-" , (41) 



as well as 




Af+l 



(42) 

X hypergeom ( [l, + 1 - 2*] , [2 + iV] , 2a;) . (43) 

This last result may be obtained by solving ([2]) recursively for xi (x) — 2x{\ ~ x) — not just by direct 
expansion of (^0)) — to obtain coefficients 

(0\k-l k-1 

^^W = ^^n(2'-^-) ■ (44) 

These A = 2 results are all well-behaved enough for the steps in the proof of the Theorem to be valid for 
X < 1/2. However, at the upper end of the interval, < x < 1/2, some additional consideration is needed. 
The expression in p3p is 1/2 at x = 1/2, independent of t and n, and therefore not small. That is to say, 
at X = 1/2 the 2~"^ prefactor in Rt (x, X,N,n) is not present to suppress the error. On the other hand, 
the ratio {dxt (x) /dx) / [dSn (x,i) /dx) in ([M)) . and in ([55t . always vanishes at x = 1/2, for any t, so the 
leading contribution to the relative error is actually zero at that point. Thus the upper end of the x interval 
does not pose a problem after all. In fact, as is evident in the graphs to follow, or by direct calculation, the 
maximum magnitude of the leading contribution to R occurs for x < 1/2, for which the 2~"^ prefactor is 
present. 

Putting all this together for the A = 2 logistic map, the leading approximation to the relative error is 

/W \ 2"-^(l-(l-2x)^"")'^^'x(l-2.x)^'(i-^"") 
i?t(x,2,iV,n)«Jw -(l-(l-2x)' ,i X ^ -. ^ . (45) 



(l-(l-2x)^') 
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For comparison purposes, we plot this last expression for various N and n, and selected t, especially to check 
(pS)) . for A = 2. That scaling law is seen to be hold fairly well, even for x « 1/2. 



0.00 

O.Oe+0 - 



-l.Oe-9 
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32^(n-5) R 




-3.0e-9 

Leading approximations to 2^("^^)i?t (x, X — 2, N — 5, n) for t = 1/2 (solid) and t ~ 3/4 (dashed) with 

n — 5, 6, & 7, in blue, red, & green, respectively. 

Finally, we note there is no discernible difference, to the accuracy of this last plot, between the leading 
approximations and the exact results for the relative error as computed from (|40p and the 5-, 6-, and 7- 
fold conjugations of '^pp''™) (^x). The largest of these differences, between the exact and the leading 
approximation of the relative error for = 5, n = 5, is shown next. 



l.Oe-14 



8.0e-15 



6.0e-15 



4.0e-15 



2.0e-15 



Delta R 



-=F- 



o.oe+o r-i — r~'~ i ~ i ~ri i i — i — i — h 

0.485 0.486 0.487 0.488 0.489 0.490 0.491 0.492 0.493 0.494 0.495 0.496 0.497 0.498 0.499 0.500 

X 

AR= Rt{x,X = 2,N = 5.n = 5)|,_, - R, (x, X = 2, N = 5,n = 5)|,,,,.„g .^^^ for t = 1/2 (solid) 

and t — 3/4: (dashed). 



So, the exact relative errors are essentially indistinguishable from their leading approximations, at least for 
A = 2 and t — 1/2 k t = 3/4, and the relative error after 7 conjugations of the 5th order series is always 
less than 3 parts in 10^^ for these two values of t. Other values of t taken from the unit interval [0, 1] are 
similarly well-approximated by the combined series and conjugation methods. 

More or less the same results can be obtained for other values of A as well, but only for one other case, 
namely A = 4, is it possible to compare to exact, closed-form results. 
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V. SUMMARY 



By combining series approximations with functional conjugations, accurate representations of continuous 
functional iterates were obtained for selected rational, sine, and logistic maps, and relative errors were esti- 
mated. These examples illustrate both the simplicity and the power of the methods. Although the rational 
and sine examples have more singular underpinnings (their Schroder functions have essential singularities) 
than the selected logistic map (whose Schroder function is analytic), nevertheless, the same methods work 
well for all three cases. 

This should be true for many other examples. Indeed, it would appear to be a relatively straightforward 
task to construct accurate approximations for the continuous iterations of all the classical special functions. 
Many of these, like the sine iterates depicted above, are expected to exhibit smooth, intuitive iterations 
that are easy to grasp, conceptually. Others, such as occur for the logistic map when the parameter leads 
to chaotic behavior, are expected to lead to much more exotic iterations. Despite the long history of and 
voluminous literature on functional equations (see P, for history and bibliographies) , the full landscape 
of features to be encountered is still largely unexplored, in our opinion. Perhaps the methods elucidated in 
this paper will help to continue that exploration. 
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